Parent Hamiltonian for the non-Abelian chiral spin liquid 



(N 

o 

(N 



o3 



i 

o 
o 



> 

m 
in 

o 

(N 



X 



Martin Greiter, 1 Darrell F. Schroeter, 2 and Ronny Thomale 3 

1 Institut fur Festskorperphysik, Postfach 3640, KIT, D 76021 Karlsruhe, Germany 
2 Department of Physics, Reed College, Portland, OR 97202, USA 
'Department of Physics, Stanford University, Stanford, CA 94305, USA 

We construct a parent Hamiltonian for the family of non-Abelian chiral spin liquids proposed re- 
cently by two of us [PRL 102, 207203 (2009)], which includes the Abelian chiral spin liquid proposed 
by Kalmeyer and Laughlin, as the special case s = i. As we use a circular disk geometry with an 
open boundary, both the annihilation operators we identify and the Hamiltonians we construct from 
these, are exact only in the thermodynamic limit. 
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Introduction. — The field of two-dimensional quantum 
spin liquids [1-15] is witnessing a renaissance of interest 
in present days [16-20]. For one thing, due to advances 
in the computer facilities available, evidence for spin liq- 
uid states in a range of models is accumulating [21, 22]. 
At the same time, spin liquids constitute the most intri- 
cate, and in general probably least understood, examples 
of topological phases [23-27], which themselves estab- 
lish another vividly studied branch of condensed matter 
physics [28-30]. If a complete description of the elec- 
tronic states in the two-dimensional (2D) CuO planes of 
high T c superconductors [31] ever emerges, the theory is 
likely based on a spin s = 1/2 liquid on a square lattice, 
which is stabilized through the kinetic energy of itinerant 
holon excitations [1]. 

Intimately related to the field of topological phases are 
the concepts of fractional quantization, and in particu- 
lar fractional statistics [32]. This field has experienced 
another, seemingly unrelated renaissance of interest in 
recent years, due to possible applications of states sup- 
porting excitations with non-Abelian statistics [33] to the 
rapidly evolving field of quantum computing and cryp- 
tography. The paradigm for this class is the Pfaffian 
state [34, 35], which has been proposed to describe the 
experimentally observed quantized Hall plateau at Lan- 
dau level filling fraction v = | [35]. The state sup- 
ports quasiparticle excitations which possess Majorana 
fermion states at zero energy [36] . Braiding of these half- 
vortices yields non-trivial changes in the occupations of 
the Majorana fermion states, and hence render the ex- 
changes non-commutative or non-Abelian [37, 38]. Since 
this "internal" state vector is insensitive to local per- 
turbations, it is preeminently suited for applications as 
protected qubits in quantum computation [39, 40]. Non- 
Abelian anyons are further established in other quantum 
Hall states including Read-Rczayi states [41], in the non- 
Abelian phase of the Kitaev model [8] , the Yao-Kivelson 
and Yao-Lee models [10, 18], and in the family of non- 
Abelian chiral spin liquid (NACSL) states introduced by 
two of us [13] . Very recently, non-Abelian statistics has 
been observed numerically in hard-core lattice bosons in 
a magnetic field, without reference to explicit wave func- 



tions [42]. 

In this Letter, we construct a parent Hamiltonian for 
the NACSL states [13]. These spin liquids support spinon 
excitations with SU(2) level k = 2s statistics for spin s, 
i.e., Abelian, Ising, and Fibonacci anyons for s = 1,1, 
and |, respectively. The method we employ here is dif- 
ferent from the method we used to identify a Hamilto- 
nian [43, 44] which singles out the Kalmeyer-Laughlin 
chiral spin liquid (CSL) state [2, 45] as its (modulo the 
two-fold topological degeneracy) unique ground state for 
periodic boundary conditions (PBCs). It is considerably 
simpler, applicable to the entire family of spin s NACSL 
states, but exact only in the thermodynamic (TD) limit 
even if we impose PBCs. 

Chiral spin liquid states. — The conceptually simplest 
way to construct the non-Abelian chiral spin liquid 
(NACSL) state [13] with spin s is to combine 2s identical 
copies of Abelian CSL states with spin 1, and project the 
spin on each site onto spin s, 



| =s® (2s- 1) ■ s-l> 



2 s 

The projection onto the completely symmetric represen- 
tation can be carried out conveniently using Schwinger 
bosons [7, 46]. For a circular droplet with open boundary 
conditions occupying TV sites on a triangular or square 
lattice, the Abelian CSL state takes the form 

l^o L ) = E ^ KL (*i, • • • , zm) S+-...- S+ m HI . . . |) 
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is a bosonic quantum Hall state in the complex "parti- 
cle" coordinates Zi = Xi + \yi supplemented by a gauge 



factor G(zi), M 



N 



i> and W are Schwinger boson cre- 



ation operators [7, 46, 47], and the w^s are those lattice 
sites which are not occupied by any of the Zi's. In this 
notation, we can write the spin s state obtained by the 
projection as 



W) = |o) 



(3) 



The lattice may be anisotropic; we have chosen the lattice 
constants such that the area of the unit cell spanned by 
the primitive lattice vectors is set to 2n. For a triangular 
or square lattice with lattice positions given by 
na + mb, where a and b are the primitive lattice vectors 
in the complex plane and n and m are integers, the gauge 
phases are simply G(j? n , m ) = (-1)(™+ I )( m + 1 ) [45, 48]. 
The NACSL state can alternatively be written as 

W> = E ^o(*i. • • • . *sn) ■ ■ ■ ■ ■ §+„ \-s) N , 

{zi,...,z„n} 

(4) 

where | — s) N = ®a=i \ s >~ s ) a i s t nc "vacuum" state in 
which all the spins are maximally polarized in the nega- 
tive z-direction, and S + are re-normalized spin flip oper- 
ators which satisfy 
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In a basis in which S z is diagonal, we may write 

1 



s 



(6) 



s - S z + 1 
Note that (5) implies 

S-(S+r\s,-s) =n(S + ) n - 1 \s,-s). (7) 

The wave function for the spin s state (3) are then ef- 
fectively given by bosonic Read-Rezayi states [41] for 
renormalized spin flips, 
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which we understand to be completely symmetrized over 
the "particle" coordinates Zj. For s = 1, they take the 
form of a Moore-Read state [34, 35] 



/ 1 \ N sN 
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For the considerations below, it is convenient to write 
the state in the form 



2 s 



m) = E #f L (*i- ■•■,**/) s+ ■...■~si M io). 
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(10) 

Since the Abelian KL CSL |V'o L ) i s an exact spin singlet 
in the TD limit TV — > oo, and is an approximate sin- 
glet for finite N, the same holds for the NACSL |i/jq) as 
well. This follows from the construction of the Schwinger 
boson projection (3), but can also be verified directly us- 
ing Perelomov's identity (see (29) in the supplementary 
material) [49] . The Abelian and non- Abelian CSL states 
trivially violate parity (P) and and time reversal (T) sym- 
metry. 

Ground state annihilation operators. — In the TD limit 
N — > oo, the NACSL ground states are annihilated by 



N 
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(S~) 2s Sp, Sl s a \r ) = 0Va, (11) 



as we will verify now. 

Let us consider the action of (S~) 2s Sg on written 
in the form (10). Since -0q l (2:i, . . . , zm) vanishes when- 
ever two arguments zi coincide, one of the Zi 's in each of 
the 2s copies in (10) must equal r\ a \ since ?/>q l (zi, . . . , zm) 
is symmetric under interchange of the z^s and we count 
each distinct configuration in the sums over {z\, . . . , zm} 
only once, we may take z\ = rj a - Regarding the action 
of Sp on (10), we have to distinguish between configura- 
tions with n = 0, 1, 2, . . . , 2s re-normalized spin flips St 
at site /?. Since the state is symmetric under interchange 
of the 2s copies, we may assume that the n spin flips are 
present in the first n copies, and account for the restric- 
tion through ordering by a combinatorial factor. This 
yields 
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where we have used (7). This implies 
n s a \r ) = {2s)\2s 
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where we have used the Perclomov identity [49], which 
states that any infinite lattice sum of e^^^ 2 G(-qp) times 
any analytic function of rip vanishes. (Strictly speak- 
ing, Perelomov [49] only considered a square lattice. The 
identity, however, holds for any 2D lattice with a sin- 
gle site per unit cell, as we show in the supplementary 
material.) 

Parent Hamiltonian. — A Hermitian, positive semi- 
definite, and translationally invariant operator which an- 
nihilates [V'o) is given by 

N 

r^^+n* = ^ Pl (s+f s (S-) 2s s+s-, (12) 

a—l ct,/3,7 

where 

W Q fl 7 = = . (13) 

This operator is not invariant under SU(2) spin rotations, 
but rather consists of a scalar, vector, and higher tensor 
components up to order 4s + 2. Since the NACSL states 
IV'o) are spin singlets, and are annihilated by T, all these 
tensor components must annihilate the state individually. 
The scalar component of T, which we denote as {r}o, pro- 
vides us with an SU(2) spin rotationally invariant parent 
Hamiltonian. 

To obtain the projected operator {r}o, we follow the 
method described in detail in ref. [50], and summarize 
here only the most important steps. With the tensor 
content of SpS~ given by 

S+S- = ^S I3 S 7 - i(Sp x S 7 y - 1%, (14) 



where 

T° 7 = A(3<^_<^s 7 ) (15) 

is the m = component of the second order tensor, we 
only need to know the scalar, vector and 2nd order tensor 
components of (S^) 2s (S~) 2s in order to obtain the scalar 
component of T. These are given by (see Sec. 5.3.2 of [50]) 

(S+) 2s (S-) 2s = a {l + aSl + 6T° Q + higher orders} 

(16) 

where 

(2*)! 2 3 V6 5 
an = , a = : b = — 7. (17 

2s + 1' s + V 2 (s + l)(2s + 3) V ' 

The scalar component of T is hence given by 

ck,/3,7 

(18) 

With S x Sp = iSp and (see Sec. 4.5.3 of [50]) 

+ 2[(S a S e ){S a S-y) + (S a S^){S a S e )], 

(19) 
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we obtain the final parent Hamiltonian 



the wave function of the NACSL (8) takes the form 



S(S + l) 2 + S a Sfi 



+ 



Pi 



2s + 3 

(S a Sp)(SaSj) + (S a S 1 )(S a Sp) 



2(8+1) 



(20) 



(It is related to (18) via {F} Q = 2a /(2s + 3) H".) This 
Hamiltonian is approximately valid for any finite disk 
with N lattice sites, and becomes exact in the TD limit 
N -> oo, where H s = 0. Note that the S a (Sp x S 7 ) 
term explicitly breaks P and T. (It would be highly de- 
sirable to identify a parent Hamiltonian which is P and T 
invariant, such that the ground states violate these sym- 
metries spontaneously, but we have so far not succeeded 
in finding one.) 

The special case s = h. — Since S +2 = for s = |, 
T™ a = for all to, and {T° a T^} Q = 0. This simplifies 
(18) significantly, and yields the parent Hamiltonian 



^ U a p^[SpS^ -iS a (Sp x S-y)] (21) 



(It is related to (18) via {r} = 2a /3 H s =i .) In 
contrast to the earlier parent Hamiltonian proposed in 
rcf. [43, 44] (SKTG) for the Abelian KL CSL (2) with 
periodic boundary conditions, (21) is not exact for finite 
N. It is considerably simpler then the SKTG model, and, 
like (20), becomes exact in the TD limit. 

Remarks on periodic boundary conditions. — It is rather 
straightforward to formulate the model on a torus. For 
simplicity, we choose the lattice constant a real, and b 
such that the imaginary part 3(6) > 0. We implement 
PBCs in both directions by identifying the sites Zj, Zi+L, 
and Zi + Lt, where L = n\a, Lt = n T a + m T b, and 
3(t) > 0. ni and m T are positive integers such that the 
number of sites N = n\m r is even, and n T is an integer. 
We place the lattice sites at positions 



with n = 0, 1, . . . , Hi — 1 and m = 0, 1, 



m=lli,_7=(m-l)M+l 
i<j 

2 "| sJV 
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where d i i (z\t) is the odd Jacobi theta function [51], and 

Zi, m = ~^2,m, (24) 



mM 



i={m-l)M+l 



are the center-of-mass coordinates and zeros, respec- 
tively. The latter can be chosen anywhere within the 
principal region bounded by the four points h[+n\a ± 
m T b), and encode the (2s + l)-fold topological degener- 
acy of the NACSL [19]. The gauge factor in (23) is given 
by 

G{r) n>m ) = (_ 1 )m.»+m e -i7r«Mm(m T -l-m) ) ^ 

where 3?(&) is the real part of b. 

The NACSL (23) is approximately annihilated by 



n: 



N 



(s-y s s i 



(26) 



for all a, where we can choose any of the three even 
Jacobi theta functions in the numerator: (u, v)=(0, 0), 
(0, ^), or (i, 0). Note that |-0g) is not strictly periodic, 
but only quasipcriodic, due to the shift of the boundary 
phases inherent in (26). The statement f2„ |-0g) « be- 
comes exact as N — > oo. 

The NACSL (23) is hence the approximate ground 
state of (20) (and for s = | also of (21)) with (13) re- 
placed by 




(27) 



where * denotes complex conjugation. As in the case with 
open boundary conditions, the model becomes exact in 
the TD limit. 

Conclusion. — We have identified a parent Hamiltonian 
for the non- Abelian CSL states [13], which becomes ex- 
act in the TD limit. This Hamiltonian should allow us 
, (22) to study the spinon and holon excitations including the 
non- Abelian braiding properties within a concise frame- 
work. The construction also extends to the Abelian s = i 
T. Then Kalmeyer-Laughlin CSL [2, 45], where it is likewise exact 



4 



only as the number of sites N — > oo, but is considerably 
simpler that the SKTG Hamiltonian [43, 44]. 
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Note added. — After this work was completed, we be- 
came aware of a manuscript by Nielsen, Cirac, and 
Sierra [52], in which they derive the s = | Hamiltonian 
(21) using null operators in the conformal correlators of 
the SU(2) level k = 1 Wess-Zumino-Witten model. 
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Supplementary material 

In this supplement, we proof the Perelomov iden- 
tity [49] for arbitrary 2D lattices using Fourier transfor- 
mation. 

The Perelomov identity. — Consider a lattice spanned 
by f] n ,m = na + mo in the complex plane, with n and m 
integer and the area of the unit cell spanned by the 
primitive lattice vectors a and b set to 2tt, 



fi = mab)\ = 2ir 



(28) 



where 9 denotes the imaginary part. Let G(ry„ im ) 
(_l)(n+i)(m+i)_ Thcn 



E 

n.m 



P(^,m)G(77„, m )e-5l""-| 2 =0 



(29) 



for any polynomial P of rj n ,m- 

Proof. — It is sufficient to proof the identity for the gen- 
crating functional 



Y d ^ m ' z G{rin,m)e-^ 



0. 



(30) 



Since G(r] n _ m ) takes the value —1 on a lattice with twice 
the original lattice constants, we may rewrite this as 



1%,. 



2 ^ e W»* e -|U».m| = Q ( 31 ) 



Kalmeyer and Laughlin [45] observed that for the square 
lattice, the second sum in (31) can be expressed as a sum 
of the Fourier transform of the function we sum over in 
the first term. We demonstrate here that their proof can 
be extended to arbitrary lattices. 

To begin with, we define the Fourier transform in com- 
plex coordinates 



/(C) = / A/faJe**^, 



(32) 



where denotes the real part and we have used (28). 
Since the area of the unit cell of our lattice is taken to be 
27r, the reciprocal lattice is given by the original lattice 
rotated by 5 in the plane without any rescaling of the 
lattice constants. In complex coordinates, 



Cn'.m' = i(n'a + m'b), 
as this immediately implies 



(33) 



Then 



= 3? ((na + mb)(— i)(n'a + m'b)) 
= nm'^s(ab) + mn'^s(ba) 
= 2tt ■ integer. 

f((n>,m')=nJ2f(Vn,m). (34) 



Eq. (34) follows directly from 

e*^'^ =Sl-£S {2) (Vn, m -v), (35) 

n' .rn' n,m 

which is just the 2D equivalent of the (Dirac comb) iden- 
tity 



E e2mn ' x = E s ( x ~ n ) 



(36) 



n — — oo 



The r.h.s. of (36) is obviously zero if x is not an integer, 
and manifestly periodic in x with period 1. To verify the 
normalization, observe that since for any N odd, 



N for y = N ■ integer 
otherwise. 



This implies 



1 ^ ^ e 27rin'y/N _ ^ 



N 



which in the limit N — > oo is equivalent to 



N N ^ 



Substituting x — y/N yields 



. 1 oo 



dx J2 e2liin ' X = 



n —-co 



(37) 



which proves the normalization in (36). 

We proceed by evaluation of the Fourier transform of 



/(C) = J dPveWe-iW e m ^ 

where we have used the integral 
d 2 r, F(n) e-sd'/l 2 -^) 

= F(adu,) J d 2 r i e-i (M2 - ilw - ,l€]) 
= F(adu,) J d 2 r]e-^ r >- w ^- w ^ 
= airF(adu,)ei wa = air F{w) 



w=0 



w=0 
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with F(ij) — e^ lz+ i 11 ^ , a = 4, and w — 2i(. If wc now substitute n' = —n, m' = —m, and hence 

Substituting (37) into (34) we obtain *'Cn',m' = Wn,m into the r.h.s. of (38), we obtain (31). 

This completes the proof. 

J2f(Vn,m) = 2 J2 e-l^'-'l 2+ ^'-^ (38) 
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